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Abstract 

This paper discusses properties of the graphs of 2-way and 3-way transportation polytopes, 
in particular, their possible numbers of vertices and their diameters. Our main results include a 
quadratic bound on the diameter of axial 3-way transportation polytopes and a catalogue of non- 
degenerate transportation polytopes of small sizes. The catalogue disproves five conjectures about 
these polyhedra stated in the monograph by Yemelichev et al. (1984). It also allowed us to discover 
some new results. For example, we prove that the number of vertices of an m x n transportation 
polytope is a multiple of the greatest common divisor of m and n. 

1 Introduction 

This paper takes a new look at the graphs of transportation polytopes. Transportation polytopes are 
well-known objects in operations research and mathematical programming (see e.g., [Tj [2"1 [51 flu ] 124 1 
[26l [271 132"! [53] and references therein). Statisticians have also interest in them of their own (see e.g., 
[HI El ESI H7J EH [25] and references therein). 

During the 1970's and 1980's the study of the classical 2-way transportation polytopes, i.e., those 
polytopes ofmxn tables satisfying row-sum and column-sum conditions, was very active. The state of 
the art of that research was carefully summarized in the comprehensive book by Yemelichev, Kovalev, 
and Kratsov [33] and Vlach's survey [32]. On the other hand, 3-way transportation polytopes, whose 
feasible points are I x m x n arrays of non-negative numbers satisfying certain sum conditions, are 
less understood. We define the polyhedra whose points are I x m x n tables satisfying certain sum 
conditions. They come in two main varieties: 
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1. First, consider the axial 3-way transportation polytope: Let x = (xi, . . . ,xi), y = (yi, . . . ,y m ), 
and z = (z\ , . . . , z n ) be three rational vectors of lengths I, m and n, respectively, with non- 
negative entries. Let P x ,y,z be the polytope defined by 1-marginals; that is, the following l+m+n 
equations in the I x ui x ti real variables (i = 1, . . . ,l;j = 1, . . . ,m; k = 1, . . . 

< a itjt k, Vi, j, k ^2 a i,j,k = x i-, Vi ai 'i' k = y i> Vj ' XT a ' 1 ^^ = Vfc - 

j,k i,k i,j 

Observe that a necessary and sufficient condition for T X) y^ z to be non-empty is that 

i j k 

and that, consequently, T x>y>z is defined by only I + m + n — 2 independent equations. 

2. Similarly, planar 3-way transportation polytopes can be defined by specifying three matrices 
U G M mjn (Q), V £ Mi >n (Q), W £ W^ m (Q) for the line-sums resulting from fixing two of 
the indices of entries and adding over the remaining index. That is, we have the following 
Im + In + mn equations (2-marginals) in the same Imn variables a>i t j,k : 

< ati^k, Vi, j, k aijfi = U jt k, Vj, k aij jk = Vi t k, Vi, k aj Jifc = Wij , Vi, j. 

One can see that in fact only Im + ln + mn — I — to — n + 1 of the defining equations are linearly 
independent for feasible systems. 

Observe that the axial 3-way transportation polytopes generalize the classical transportation poly- 
tope of size to x Ti, which coincides with T Xi y iZ for I = 1 and x = ^2Vj = z k- A less trivial rewriting 
of the classical 2 x n transportation polytope as a 2 x 2 x n 3-way planar transportation polytope is 
given in Theorem 11.41 below. 

Recall that the 1-skeleton or graph of a convex polytope P is the set of all O-dimensional and 
1-dimensional faces (vertices and edges) of P, with their natural incidence relation. The main focus 
of this paper is to investigate the number of vertices and the diameters of the graphs of classical (that 
is, 2-way) and 3-way transportation polytopes. 

Some of the statements below require our transportation polytopes to be non-degenerate. By this 
we mean that the polytope is simple (i.e., every vertex is adjacent to dimension many other vertices) 
and it is of maximal possible dimension (that is, dimension Imn — I — m — n + 2 for I x to x n axial 
transportation polytopes, and dimension (/ — l)(m — l)(n — 1) for I x to x n planar transportation 
polytopes). Graphs of non-degenerate transportation polytopes are of particular interest because they 
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have the largest possible number of vertices and largest possible diameter among the graphs of all 
transportation polytopes of given type and parameters. Indeed, if a transportation polytope P is 
degenerate, by carefully perturbing the marginals that define it we can get a non-degenerate one P' . 
The perturbed marginals are obtained by taking a feasible point X in P, perturbing the entries in the 
table and using the recomputed sums as the new marginals for P' . The graph of P can be obtained 
from that of P' by contracting certain edges, which cannot increase either the diameter nor the number 
of vertices. 

Our main result is a bound on the diameter of axial 3-way transportation polytopes: 

Theorem 1.1 The graph of the I x m x n axial transportation polytope T xy>z has diameter at most 
2(l + m + n- 3) 2 . 

A similar result for the graph of a classical transportation polytope was given by Brightwell et al. 
(see [7j), who proved an upper bound of 8(m + n — 2) for the diameter. More recently, Hurkens (see 
|20j) has obtained a bound of 4(m + n — 1), a factor of four away from the predicted value of the 
Hirsch conjecture. 

Using computational tools, we also give a complete catalogue of non- degenerate 2-way and 3-way 
transportation polytopes (both axial and planar) of small sizes. This allowed us to explore properties 
of transportation polytopes (e.g., their diameters and how close they were to the Hirsch conjecture 
bound). The summary of the catalogue is: 

Theorem 1.2 • The only possible numbers of vertices of non- degenerate 2x3, 2x4, 2x5, 3x3, 
and 3x4 classical transportation polytopes are those given in Table [7J 

• The only possible numbers of vertices of non- degenerate 2x2x2, 2x2x3, 2x2x4, 2x2x5, 
and 2x3x3 planar transportation polytopes are those given in Tabled 

Every non- degenerate 2x3x4 planar transportation polytope has between 7 and 480 vertices. 

• The only possible numbers of vertices of non- degenerate 2x2x2 and 2x2x3 axial transportation 
polytopes are those given in Tabled 

Every non- degenerate 2x2x4 axial transportation polytope has between 32 and 504 vertices. 
Every non- degenerate 2x3x3 axial transportation polytopes has between 81 and 1056 vertices. 
The number of vertices of non- degenerate 3x3x3 axial transportation polytopes is at least 729. 

The catalogue was obtained via the exhaustive and systematic computer enumeration of all com- 
binatorial types of non-degenerate transportation polytopes. The theoretical foundations of it are the 
notions of parametric linear programming, chamber complex, Gale diagrams and secondary polytopes 
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Size 


Dimension 


Possible numbers of vertices 


2x3 


2 


3 4 5 6 


2x4 


3 


4 6 8 10 12 


2x5 


4 


5 8 11 12 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 


3x3 


4 


9 12 15 18 


3x4 


6 


16 21 24 26 27 29 31 32 34 36 37 39 40 41 
42 44 45 46 48 49 50 52 53 54 56 57 58 60 61 
62 63 64 66 67 68 70 71 72 74 75 76 78 80 84 90 96 



Table 1: Numbers of vertices possible in non-degenerate classical transportation polytopes 



Size 


Dimension 


Possible numbers of vertices 


2x2x2 


1 


2 


2x2x3 


2 


3 4 5 6 


2x2x4 


3 


4 6 8 10 12 


2x2x5 


4 


5 8 11 12 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 


2x3x3 


4 


5 8 9 11 12 13 14 15 16 17 18 19 20 21 22 23 
24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 
41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 



Table 2: Numbers of vertices possible in non-degenerate planar transportation polytopes 



Size 


Dimension 


Possible numbers of vertices 


2x2x2 


4 


8 11 14 


2x2x3 


7 


18 24 30 32 36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 
68 70 72 74 76 78 80 84 86 96 108 



Table 3: Numbers of vertices possible in non-degenerate axial transportation polytopes 
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(see |15j). We present them in Section [2j The full catalogue of transportation polytopes (includ- 
ing other families, such as 3 x 5, 4 x 4, 4 x 5, etc.) is available in a searchable web database at: 
http : // www. math. ucdavis . edu/~ekim/transport at ion_polytope_dat abase/. 

Based on the data we collected, we discovered and proved (in Section [3} the following results: 

Theorem 1.3 The number of vertices of a non- degenerate m x n classical transportation polytope is 
divisible by GCD(m,n). 

Theorem 1.4 The 2 x 2 x n planar transportation polytopes are in 1 — 1 correspondence with the 2 x n 
classical transportation polytopes, with corresponding pairs being linearly isomorphic. 

Note that Theorem 11.41 is best possible in the sense that for m, n > 3 there are many more types 
of planar 2 x m x n transportation polytopes than types ofmxn transportation polytopes - see the 
rows of the 3x3 and 2x3x3 polytopes in Tables [1] and [2j Finally, we state two open conjectures at 
the end of Section [2j 

We close the introduction explaining the relevance of our results to current research. Bounding 
the diameter of graphs of polytopes has received a lot of attention because of its connection to the 
performance of the simplex method for linear programming and, most especially, to try to understand 
the Hirsch conjecture (see |23|. 122] and references therein). 

We recall that the Hirsch conjecture asserts that the diameter of the graph of any polytope of 
dimension d and with / facets is bounded above by / — d. Not only is this conjecture open, but even 
the weaker statement asserting a polynomial upper bound (in / and d) for the diameters of graphs of 
polytopes is unknown (although a quasi-polynomial bound appeared in |22j). 

As we observed, [7] provided the first linear bound for the diameter of the graphs of 2-way trans- 
portation polytopes. Theorem 11.11 provides a quadratic one for axial 3- way transportation polytopes 
and, moreover, a sublinear one if we assume that the three parameters /, m and n are approximately 
the same. (Observe that the number of facets of a 3-way transportation polytope is bounded above 
by the product Imn of its size parameters). 

Bounding the diameter of 3-way transportation polytopes is particularly interesting because of the 
following results recently proved by two of the authors in |13j : 

1. Any rational convex polytope can be rewritten as a face F of an axial 3-way transportation 
polytope. The sizes I, m, n, 1-marginals x, y, z, and the entries aij t k that are prescribed to be 
zero in the face F can be computed in polynomial time on the size of the input. 

2. More dramatically, any convex rational polytope is isomorphically representable as a planar 
3-way transportation polytope. 
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That is to say, a version of Theorem 11.11 for the 3- way planar case, or a version for the axial case 
that allows one to prescribe some variables to be zero, would provide a polynomial upper bound on 
the diameter of the graph of every convex rational polytope. 

Another consequence of these results is that the method of Section [2] for enumerating all combi- 
natorial types of planar 3-way transportation polytopes, yields, in particular, an enumeration of all 
types of rational convex polytopes. 

Let us finally mention that our systematic listing of non-degenerate transportation polytopes pro- 
vides the solution to at least four open problems and conjectures about transportation polytopes 
stated in the monograph [33j : 

1. Klee and Witzgall in [24] prove that the largest possible number of vertices in classical transporta- 
tion polytopes of size mxnis achieved by the generalized Birkhoff polytope (the transportation 
polytope with parameters X{ = n V i, yj = m V j). Problem 32 in page 400 of [33] conjectured 
that the same holds in general. 

But in Example 12. 71 we provide an explicit counterexample of this for planar 3-way transportation 
polytopes. (In this case, the generalized Birkhoff polytope is the planar 3-way transportation 
polytope whose 2-marginals are given by the mxn matrix U(j, k) = I, the I x n matrix V(i, k) = 
m, and the I x m matrix W(i,j) = n.) 

2. Question 36 on page 396 of [33] asked: Is it true that every integer of the form (I — l)(m — 
l)(n — 1) + 1 where 1 < t < ml + nl + ran — I — m — n, and only these integers, can equal the 
number of facets of a non- degenerate planar 3-way transportation polytope of order I x m x n, 
where l,m,n > 2? 

For the case I = m = 2 and n = 3, the conjecture asks if every integer from 3 to 11, and only 
these integers, equal the number of facets of non-degenerate 2x2x3 planar transportation 
polytopes. Since in this case the number of facets equals the number of vertices (because the 
polytopes are two dimensional) Table answers the question negatively: only facet-counts from 
3 to 6 occur, while 7 through 11 are in fact missing. 

3. Similarly, Conjecture 33 on page 400 of [33J asked: Is it true that every integer from 1 to 
ml + nl + mn — I — m — n + 1, and only these numbers, are realized as the diameter of a planar 
3-way transportation polytope of order I x m x n? 

The same case I = m = 2, and n = 3 shows that this is false. The transportation polytopes 
obtained are polygons with up to six sides, hence of diameter at most three, instead of 10. 
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4. Open problem 37 in page 396 of [33J asks whether the numbers of vertices of I x m x n non- 
degenerate planar transportation polytopes satisfy: 

(I - l)(m - l)(n - 1) + 1< f < 2(1 - l)(m - l)(n - 1). 

We show the answer is no even in the case 2x2x4. 

In addition to the four solved problems above, Theorems 11.21 and 11.31 are initial steps on the 
solution of Problem 25 in page 399 of [33] . It asks to find the complete distribution of possible number 
of vertices for transportation polytopes. 

2 Classifying Transportation Polytopes 

Theorem 11.21 was obtained through an exhaustive enumeration whose foundation is the theory of 
secondary polytopes and parametric linear programming. In some cases when the full enumeration 
was impossible we at least get lower and upper bounds for the number of vertices that these polytopes 
can have. In this section we discuss the necessary background to understand the construction of the 
complete catalogue. 

2.1 Enumeration via Regular triangulations and Secondary polytopes 

We begin by recalling some basic facts about convex polytopes presented, as all transportation poly- 
topes are, in the form P c = {x : Bx = c,x > 0}. For the case of transportation polytopes, the vector 
c is the vector given by the demand/supply quantities. Fix a matrix B of full row rank. Most of our 
results are obtained by studying what happens to the combinatorics of P c as the vector c changes while 
we fix the matrix B. This study, for general matrices, is known as parametric linear programming (see 
Chapter 1 of [IS]). 

A subset of W 1 that is closed under addition and under multiplication by positive scalars is a cone. 
For any set L of vectors in W 1 , the cone generated by L, denoted cone(L), is the set of all vectors that 
can be expressed as non-negative linear combinations of the members of L. Abusing notation, for a 
matrix B, by cone(-B) we mean the cone generated the set of column vectors of B. 

A maximal linearly independent subset b of B is a basis of B. Geometrically, each basis of the 
matrix B spans a simple cone inside cone(-B). Every basis b of B defines a basic solution of the system 
as the unique solution of the m linearly independent equations bxb = c and Xj = for j not in b. A 
basic solution is feasible if in addition, x > 0. Geometrically, a basic feasible solution corresponds to 
a simple cone that contains c. In fact, one can see that P c is non-empty if and only if c € cone(-B). 
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A fundamental fact in linear programming is that, for a given right-hand-side vector c, all vertices 
of the polyhedron P c are basic feasible solutions (see [301 S3]). Moreover, if the polyhedron P c is 
assumed to be non-degenerate then the basic feasible solution must be strictly positive on the entries 
corresponding to the basis b. Geometrically, a basis d C B produces a vertex of a non-degenerate 
P c if and only if c lies in the interior of the cone generated by b. In conclusion, the vertices of a 
non-degenerate polyhedron P c , from this family of parametric polytopes, are in bijection with the 
bases that contain the right-hand-side vector c within their interior. 

We now look at what happens when we let c vary. If P c is non-degenerate and the change in c is 
small, the facets of P c move but the combinatorial type of P c does not change. Only when a basic 
solution changes from being feasible to not feasible, or viceversa, the combinatorics of P c (that is, the 
face lattice and, in particular, the graph of P c ) can change. 

Put differently: Denote by the set of all cones generated by bases of B. Let <9X_b denote 
the union of the boundaries of all elements of The connected components of cone(-B) \ 8T,b are 
open convex cones called the chambers of B. Equivalently, we call the chamber associated to a given 
right-hand side vector c the intersection of the interiors of simple cones that contain c in their interior. 
We remark that every feasible and sufficiently generic c is in a chamber (as opposed to lying on dY,^). 
Two vectors c\ and C2 in the same chamber determine non-degenerate polytopes P Cl and P C2 that are 
equivalent up to combinatorial type. The collection of all the chambers is the chamber complex or 
chamber system associated with B. Putting all this together we conclude 



Proposition 2.1 To represent all the possible combinatorial types of non-degenerate polytopes of the 
form P c = {x : Bx = c, x > 0}, for a fixed B and varying vector c, it is enough to choose one c from 
each chamber of the chamber complex of B. 



Example 2.2 Consider the matrix i?2,3; the constraint matrix of all 2 x 3 transportation polytopes. 
That is: 

r 1 1 1 



B 



2,3 



1 1 1 

10 10 

10 10 

1 1 



This means, the system {-B2.3 V = c, y > 0} defines the 2x3 transportation polytopes with marginals 
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The columns of the matrix #2,3 span a four-dimensional cone in Mr. It will be relevant later 
that if we slice this cone by an affine hyperplane (such us YlVi = 1) we obtain the three-dimensional 
triangular prism shown in Figure [7J but embedded in M 4 . 

The chamber complex can be obtained by slicing the prism with the six planes containing a vertex 
of the prism and the edge "opposite " to it. The resulting chamber complex is hard to visualize or draw, 
even in this small case, but we will see below how to recover the structure of the chamber complex 
for this example using Gale transforms. In particular, as we will see, this decomposes the triangular 
prism into 18 chambers. 

It is very easy to "sample" inside the chamber complex and find chambers of different numbers of 
bases, i.e., transportation polytopes with different number of vertices. One can simply throw random 
positive values to the cell entries of an I x m x n table and then compute the 1-marginals or 2-marginals 
associated to it. But with this method it is not obvious how to guarantee that one has obtained all the 
possible chambers. For this we use the approach based on Gale transforms and regular triangulations, 
that we now explain. 

A vector configuration A of r vectors in R( r - d ) is called a Gale transform of another vector 
configuration B of r vectors in M. d if the row space of the matrix with columns given by A is the 
orthogonal complement in W of the row space of the matrix with columns given by B. Gale transforms 
are essential tools in the study of convex polytopes because the combinatorial properties of B and A 
are intimately related (see Chapter 6 in |34j for details). Proofs of the following statements can be 
found in [5J [19] . See also Chapters 4 and 5 in [15] . 

Lemma 2.3 • The chambers of B are in bisection with the regular triangulations of the Gale 
transform B of B: From a chamber in B, one can recover a regular triangulation of B via 
complementation, namely for a basis a of vectors in B the elements of B not belonging to a 
form a basis for B. The collection of those bases gives a triangulation of B (see below for an 
example). 

• There exists a polyhedron, the secondary polyhedron, whose vertices are in bijection with the 
regular triangulations of the Gale transform B. 

• The face lattice of the chamber complex of the vector configuration B is anti-isomorphic to the 
face lattice of the secondary polyhedron of the Gale transform B of B. The latter is, in turn, 
isomorphic to the refinement poset of all regular subdivisions of B. 

• If B defines a pointed polyhedral cone (for example, if all its entries are non-negative as it 
is the case for transportation polytopes), then its Gale transform B is a totally cyclic vector 
configuration. That is, the cone spanned by B is the whole of W 1 . 
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Example 2.4 (Example 12.21 continued) We take again the matrix i?2,3 (the case of 2 x 3 trans- 
portation polytopes). A Gale transform consists of the columns of the matrix 

1-1 -110 
1 -1-10 1 

In Figure{l\we represent the Gale diagram £?2,3 and its 18 regular triangulations, each one providing 
a combinatorial type of non- degenerate 2x3 transportation polytope, although repeated combinatorial 
types occur. The chamber adjacency, which corresponds to bistellar flips, is indicated by dotted edges. 

1 3 



B 



2,3 




Figure 1: The 18 regular triangulations of i?2,3i which correspond to the 18 combinatorial types of 
non-degenerate 2x3 transportation polytopes. 



Thus, generating all the combinatorial types of non-degenerate transportation polytopes is the 
same as listing the distinct regular triangulations of the Gale transform of the defining matrix B. Note 
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that in our case B depends only on the sizes l,m,n and the type, axial or planar, of transportation 
polytopes we look at. That is, we have one B for each row of Tables [H [2] and [3l 

Now, it is well-known that the regular triangulations of a vector configuration can all be generated 
by applying bistellar flips to a seed regular triangulation (see [U [I5j [34]). Bistellar flips are combina- 
torial operations that transform one triangulation into another and regularity of triangulations can be 
determined by checking feasibility of a certain linear program (in our case, the very one that defines 
P c ). An example of the linear programming feasibility problem is given in Chapter 5 of |15j . 

Example 2.5 (Example 12.21 continued) Consider the only triangulation of B2 3 with six cones (the 
first one of the middle row in Figure^ . The necessary and sufficient conditions in the non-negative 
vector (ci,C2, ■ ■ ■ ,cq) in order to produce this triangulation are that each Ci be smaller than the sum of 
the two adjacent to it. That is, 

ci<c 5 + c 6 , c 2 <c 4 + c 6 , c 3 <c 4 + c 5 , 

C4<C2+C3, C5<Cl+C3, Cq < C\ + C2- 

Thus, these conditions on the marginals characterize the 2x3 transportation polytopes that are 
hexagons. 

To implement this method we have written a C ++ program that is available from the web page of 
the second author. This program calls TOPCOM (see |29j). a package for triangulations that computes, 
among other things, the list of all regular triangulations of a configuration. Our program also calls 
polymake (see |18j ) for the Gale transform, and computes one vector c per chamber. The output is a 
list of transportation polytopes, one per chamber, given in the polymake file format. 

2.2 Lower and upper bounds via integer programming 

Even for seemingly small cases, such as 3 x 3 x 3 transportation polytopes, listing all chambers (and 
thus all combinatorial types of transportation polytopes) is practically impossible. In these cases 
we have followed a different approach to at least obtain upper and lower bounds for the number of 
vertices of transportation polytopes. By the discussion above, this is the same as finding bounds 
for the number of simplices in triangulations of the Gale transform B. Here we follow the method 
proposed in [11], based on the universal polytope. This universal polytope, introduced by Billera, 
Filliman and Sturmfels in [3], has all triangulations (regular or not) of a given vector configuration A 
in W 1 as vertices, and projects to the secondary polytope. The universal polytope has much higher 
dimension than the secondary polytope, in fact its ambience dimension is the number of possible bases 
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of the configuration A, thus no more than (J A) • It has the advantage that the number of simplices 
in different triangulations is given by the values of a linear functional ip. 

More precisely, we think of the chambers of cone(B) as the vertices of the following high-dimensional 
0/1-polytope: Assume B is a vector configuration with n vectors inside M. d . Let N be the number of 
(i-dimensional simple cones in B. We define Ub as the convex hull in ~R. N of the set of incidence 0/1 
vectors of all chambers of B. For a chamber T the incidence vector vt has coordinates (ut)o- = 1 if 
the basis a G T and (t>r)o- = if a is not a basis of T. The polytope Ub is the universal polytope 
defined in general by Billera, Filliman and Sturmfels in [3] (although there it is defined in terms of 
the triangulations of the Gale transform of B). 

In [11], it was shown that the vertices of the universal polytope of B are exactly the integral points 
inside a polyhedron that has a simple inequality description in terms of the oriented matroid of B 
(see |11| [34] for information on oriented matroids). The concrete integer programming problems in 
question were solved using C-plex Linear Solver™ . The program to generate the linear constraints 
is a small C ++ available from the web page of the first author (see |31j). 

Example 2.6 (Example 12.21 continued) Continuing with the running example, if B is i?2,3 from 
above, then Ub is defined in M 15 , where each coordinate is indexed by a 2-subset a of {1, . . . ,6}. Let 
S denote the set of all bases. Then N = \S\ = 15. Thus Ub is the convex hull in M. N of the incidence 
vectors corresponding to the 18 chambers of B. By Lemma \2.c\ this is equivalent to the convex 
hull of the incidence vectors vt of the 18 triangulations of B. For example, the triangulation T = 
{{1, 2}, {1, 3}, {2, 4}, {3, 4}} in Figure [7] gives the incidence vector vt = ^{1,2} + e {i,3} + e {2,4} + e {3,4} 
(where e a is the basis unit vector in the direction a) as one of the vectors of the convex hull. 

The convex hull of these 18 incidence vectors is a 6 -dimensional 0/1 polytope Ub in M 15 . That the 
dimension is (at most) six follows from the following considerations: 

• Since the pairs {1,4}, {2,5} and {3,6} are not full- dimensional and thus never appear as a 
simplex in any triangulation T of B, Ub is contained in the subspace ^{1,4} = £{2,5} = £{3,6} = 0- 

• Since the vector 1 has 2 and 6 on one side and 5 and 3 on the other, in every triangulation the 
sum ^{1,2} + £{1,6} equals the sum £{1,3} + £{1,5} (and it equals zero or one depending on whether 
the triangulation uses the vector 1 or not). This implies the first of the following equalities, the 
rest being the analogue statement for the other five vectors. 

£{1,2} + £{1,6} ~ £{1,3} - #{1,5} = 
£{2,3} + ^{2,4} - £{1,2} - £{2,6} = 
£{1,3} + £{3,5} - £{2,3} - £{3,4} = 
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x {3,4} + £{4,5} ~ £{2,4} ~ x {4,6} = 
£{1,5} + £{5,6} ~ £{3,5} ~ £{4,5} = 
£{2,6} + £{4,6} ~ £{1,6} ~ x {5,6} = 

Observe that one of these equations is redundant, since the sum of the left-hand sides is already 
zero. 

• Since every triangulation needs to cover the angle between, for example, vectors 1 and 6, and 
this angle is covered only by the cones 16, 12 and 56 (see FigureUty, we have that 

^{1,6} + ^{1,2} + £{5,6} = !• 

The results in \1V§ say that Ub is the convex hull of the non-negative integer points in M 15 satisfying 
this list of equations. 

We now denote by ifi € (M N ^j * the cost vector defining the linear function 

ijj(x) = (1,1,... ,1) ■ x = J2 X ° 

Then the values of ip(x) on Ub H {0,1}^ are the only possible values for the number /o of vertices 
of non- degenerate polytopes of the form P c = {x\Bx = c, x > 0}. In particular, the solutions to 
the linear programming relaxations: "minimize (respectively maximize) ip(x) subject to x € Ub" give 
lower (respectively upper) bounds to the possible values for the number /o of vertices of non- degenerate 
polytopes of the form P c = {x\Bx = c, x > 0}. In the running example, 3 < ip(x) < 6 whenever 
x € Ub H {0, 1}^. From Tabled we observe that the number f of vertices of a non- degenerate 2x3 
transportation polytope equals 3 ; 4, 5 or 6. 

Example 2.7 Here is an application of our method. Table^is an explicit vector of 2-marginals for a 
3x3x3 transportation polytope which has more vertices (270 vertices) than the generalized Birkhoff 
polytope, with only 66 vertices. 



164424 


324745 


127239 


163445 


49395 


403568 


184032 


123585 


269245 


262784 


601074 


9369116 


1151824 


767866 


8313284 


886393 


6722333 


935582 


149654 


7618489 


1736281 


1609500 


6331023 


1563901 


1854344 


302366 


9075926 



Table 4: Counterexample to [33, open problem 37] 



Based on the data collected from the enumeration process, we also conjecture to be true: 
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Conjecture 2.8 The graph of every non- degenerate mxn transportation polytope has a Hamiltonian 
cycle (mn > 4). 

Conjecture 2.9 If P is a non- degenerate I x m x n axial transportation polytope (l,m,n >3), then 
the diameter of its graph G(P) is equal to f — d, where d = Imn — I — m — n + 2 is the dimension and 
f is the number of facets of P. 

For mxn classical transportation polytopes (m, n < 5), there are non-degenerate polytopes where 
the diameter of the graph G(P) is strictly less than f — d. 



3 Proofs of Theorems 11.31 and 11.41 

We start with Theorem 1 1.41 The 2 x 2 x n planar transportation polytopes are linearly isomorphic to 
the 2 x n classical transportation polytopes. 

Lemma 3.1 The planar 2 x m x n transportation polytopes are exactly the mxn transportation 
polytopes with bounded entries. 

Proof. Every planar 2 x m x n transportation polytope 

P = \ i a i,j,k) G K^ xmxn : dij t k = 0>i,j,k = Ui,k, aij,fc + 02,j,fc = 

{ k j 

is linearly isomorphic to an m x n transportation polytope with bounded entries, 

Q = { ( a l,j,k) ^ K+ Xn : a l,j,k = x l,j, 5Z a l,j,A; = Ul,k, a l,j,k < Zj t k 

(hi 

via the projection K 2xmxn — > K mx " taking (a,i j t k) i— > (aij^), which maps P bijectively onto Q. 
Conversely, every mxn transportation polytope Q with bounded entries is linearly isomorphic to a 
planar 2 x m x n transportation polytope P by defining X2j '■= (Ylk z j,k) ~ x i,i> 3 = 1> • • • > m an( ^ 
V2,k ■= (Ej z j,k) ~ Vl,k, k = l,...,n. ZD 



Proof of Theorem \l-4\ Consider any planar 2 x 2 x n transportation polytope 

, a i,j,k = a i,l,k + 0>i,2,k = Ui,k, a l,j,k + a 2,j,k = z j,k 



I k 



14 



The equations of the last two types imply that for each k we can express all the a,ij k in terms of a\ \ k 
as follows: 



°i,2,fc — Vl,k — a l,l,Jfc) 
02,i,fc = Zi.fc — ai,i,fc, 



In particular, P is linearly isomorphic to its projection 

Q = |(ai,i, fc ) £ M^ x2x ™ : a fc < a 1Xk < /3*,J^ai, llfc = 2:1,1 1 , 

where a fc = max{0,zi ifc - y 2yk \ = max{0,yi jfc - z 2:k } and (3 k = min{yi ifc , z 1>k }. Now, by applying a 
translation to Q, there is no loss of generality in assuming that a k = for all k. Then Q is a 1-way 
transportation polytope with bounded entries, isomorphic (by Lemma l3.1j) to a 2 x n transportation 
polytope. 

Conversely, any 2 x n transportation polytope 



Q = (aj, k ) G M^. xn 



: ^ a J>fc = Xj, a Xife + a 2ifc = y fc > , 
A; J 



is linearly isomorphic to the following planar 2 x 2 x n transportation polytope: 



: 2 L'xn . Sfe a l,j,k — J2k a 2,3-i,fc — Xj, 



Z)i a «,l,fc = Si a «,2,fc = a l,i,fc = Y^j a 2,j,k = Uk J 

The equations relating the solutions of Q to those of P are aj tk = a>i t j tk = a-2,s-j,k- H 

One final comment. The above result is best possible since the list of 2 x 3 x 3 planar transportation 
polytopes presented in Table [2] is not the same as the list of 3 x 3 classical transportation problems 
presented in Table [TJ 

We now move to Theorem 11.31 The number of vertices of a non- degenerate m x n classical trans- 
portation polytope is divisible by GCD(m,n). The first observation, already hinted in Example 12.21 is 
that the vector configuration B m ^ n associated to these transportation polytopes is (a cone over) the set 
of vertices of the product A mjT1 of two simplices of dimensions m — 1 and n — 1. So, we are interested 
in the cardinalities of chambers in the product of two simplices. Here and in what follows we call 
the cardinality of a chamber c of B m ^ n the number of bases of B m ^ n that contain the chamber c. We 
denote it by \c\. The proof of Theorem 11.31 consists of the following two steps, which are established 
respectively in the two lemmas below: 
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• There is a "seed" chamber in A m ^ n whose cardinality is indeed a multiple of GCD(m, n). 

• The difference in the cardinalities of any two adjacent chambers of A m x A n is a multiple of 
GCD(m,n). 

Since the chamber complex is a connected polyhedral complex (where two adjacent chambers are 
divided by a hyperplane supported on the vector configuration) the two lemmas settle the proof. 

Let us define the lexicographic chamber of A m) „ recursively as the (unique) chamber incident to 
the lexicographic chamber of A mj „_i. The recursion starts with A mj i, which is an (m — l)-simplex and 
contains a unique chamber. Observe that the definition of the lexicographic chamber is not symmetric 
in m and n. For example, the lexicographic chamber of the triangular prism A^^ is the one incident 
to a basis of the prism, and has cardinality 3. The lexicographic chamber of A2,3 is incident to one of 
the edges parallel to the axis of the prism, and has cardinality four. 

Lemma 3.2 The lexicographic chamber of A m ^ n is contained in exactly m n ~ l simplices. 

Proof. The cardinality of the lexicographic chamber of A m ^ n equals the cardinality of the lexicographic 
chamber of A mjn _i times the number of vertices of A mj „ not lying in its facet A mj „_i. The latter 
equals m. H 



When moving from a chamber c_ to an adjacent one c+ we "cross" a certain hyperplane H spanned 
by all except one of the elements of any basis containing c+ but not containing c_. Let us denote by 
(7+ and C_ the subsets of B m ^ n lying in the sides of 7i containing c+ and c_ respectively (Remember 
that, in our case, B m<n equals the set of vertices of A mjn ). Observe also that the common boundary 




Figure 2: A cross-section of the chamber complex of some cone with two adjacent chambers 



Co C H of c + and c_ is a chamber in the vector configuration B m ^ n fl H. 
Lemma 3.3 1. |c+| - |c_| = |co|(|C+| - |C_|). 
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2. If B mn is the set of vertices of A mtU , then \C+\ — |C_| is a multiple of GCD(m,n). 

Proof. A basis b + contains c + but not c_ if and only if b + is of the form 6q U where 60 is a basis 

of B n TL containing cq and v+ is an element of C+. This and the analogous property for c_ proves 
the first part. 

For the second part, we restate a few facts in the terminology of oriented matroids. This makes 
the proof easier to write (for details we recommend [6] ) : 

• In oriented matroid terminology a pair (C+, C_) consisting of the subconfigurations on one and 
the other side of a hyperplane H spanned by a subset of B m>n is called a cocircuit of B m ^ n . That 
is, part 2 is a statement about the cocircuits in the oriented matroid hAm n 

associated to the 

vertices of the product of two simplices. 

• The oriented matroid Ai mn coincides with the one associated to the complete directed bipartite 
graph K mn . (i.e., the complete bipartite graph with all of its edges oriented from one part to the 
other). Thus, part 2 is a statement about the cocircuits in the oriented matroid of the directed 
K 

• The cocircuits of a directed graph G = (V, E) are all read off from cuts in the graph. By this 
we mean that the vertex set V is decomposed into two parts (V+, VI). The cocircuit (C+, C_) 
associated to the cut (V+, V_) has C+ consisting of all the edges directed from V+ to V- and C_ 
consisting of all the edges directed from V_ to V+. 

Using the dictionary between the directed graph K m ^ n and the product of simplices we can finish 
the proof. Let (V + , V_) be a cut in the complete directed bipartite graph K m n . Since our graph is 
bipartite, we have V+ and V_ naturally decomposed as V+. U and U respectively. The 
sizes of C+ and C_ are then: 

|C+| = |vJ m) |-|vi n) | and |C_| = |vi m) |-|v| n) |. 

Now, using that | v| m ^ | + \V_ \ = m and |vl | + |V2 | = n we get: 

t /~1 I \ lT/( m )l I lT^( n )l\ lT/( n )l I lT^( m )l\ lT^( m )l lTA( n )l 

|C+| — I C — I = |V_^ I • (n — I Vj: I) — I V+ I • (m — \V+ \) = \Vj r \ ■ n — \V+ J \ ■ m, 
which is clearly a multiple of GCD(m, n). H 
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4 The diameter of 3-way axial transportation polytopes 



Here we consider a 3-way axial transportation polytope T x ,y,z denned by certain 1-marginal vectors 
x, y and z. Recall that for bounding its diameter there is no loss of generality in assuming T X)VyZ 
non-degenerate, that is, that x, y and z are sufficiently generic. In the non-degenerate case, at every 
vertex V of our polytope exactly Imn — I — m — n + 2 variables are zero, and exactly I + m + n — 2 
are non-zero. The set of triplets (i,j,k) indexing non-zero variables will be called the support of the 
vertex V. 

We say that a vertex V of T xyz is well-ordered if the triplets (i,j,k) that form its support are 
totally ordered with respect to the following coordinate- wise partial order: 

k) < (i',f, k') if i < i' , and j < f, and k < k'. (1) 

Observe that a set of I + m + n — 2 triplets satisfying this must contain exactly one triplet k) 
with i + j + k = p for each p = 3, . . . ,1 + m + n. Actually, supports of well-ordered vertices are the 
monotone staircases from (1, 1, 1) to (I, m, n) in the I x m x n grid. 

Lemma 4.1 If x, y and z are generic, then T x ,y,z has a unique well-ordered vertex V. 

Proof. Existence is guaranteed by the "northwest corner rule algorithm" , which fills the entries of the 
table in the prescribed order (see survey [28J or exercise 17 in Chapter six of [33]). More explicitly: let 
Vi,m,n = min{x;, y m , z n }. Genericity implies that the three values xi, y m and z n are different. Without 
loss of generality we assume that the minimum is z n . Then, our choice of Vi^ m ,n makes Vij iTl = for 
every other pair (i,j). The rest of our vertex V is a vertex of the (xmx(n-l) axial transportation 
polytope with margins x' = (x\,. . . xi~z n ), y' = (yi,.. . , y m -i, Vm-Zn), and z' = (zi, . . . ,z n -x). 

Uniqueness follows from the same argument, simply noticing that the support of a well-ordered 
vertex always contains the entry (l,m,n), and no other entry from one of the three planes (/,*,*), 
(*, m, *) and (*, *, n). This, recursively, implies that the vertex can be obtained by the northwest 
corner rule. H 

Remark 4.2 Another proof of Lemma \4-1\ can be done using the formalism of chambers developed in 
the previous sections: it is obvious (and is proved in Ulf ) that if c denotes a chamber of B, and T 
is a triangulation of cone(B), then there is a unique maximal- dimension simplex in T that contains 
c. Thus, Lemma follows from the fact that monotone staircases in the I x m x n grid form a 
triangulation of the vector configuration Bi jm>n of axial I x mx n transportation polytopes. The latter 
is well-known, once we observe that Bi mn is the vertex set of a product of three simplices. The 
triangulation in question is called the "staircase triangulation" of it (see Chapter 6 of 
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Example 4.3 To illustrate Lemma consider the non- degenerate 3x3x3 axial transportation 
polytope T X ^ )Z with: 



y~] ai,j,k = 112 
Qi,l,fe = 40 

i,k 

^ ^ a, j i = 82 

i,3 



y~] a 2,j,fc = is 

3,k 



^ a3,i,fe = 30 

^ Oi,2,*i = 6 ^ ai )3) fc = 114 
^ ai ,i,2 = 44 ^ a itjt3 = 34 



The ufiiQUG well- ordered vertex V of T x ^y^ z has the non-zero coordinates ct(i,i,i) = 40, a(i,2,i) = 6, 
a (i,3,i) = 36, 0(1,3,2) = 30, 0(2,3,2) = 14, a (2,3,3) = 4, and 0(3,3,3) = 30. Note that the non-zero entries 
of V are totally ordered (they are presented above in increasing order) with respect to Figure^ 
depicts the associated monotone staircase. 




Figure 3: A well-ordered vertex and its staircase. 



Our bound on the diameter of T x> y jZ is based on an explicit path that goes from any initial vertex 
V of T XjV)Z to the unique well-ordered vertex V. To build this path we rely on the following stratified 
version of the concept of well-ordered vertex. We say that a vertex V of T XtV)Z is well-ordered starting 
at level p, where p is an integer between 3 and I + m + n if: 

1. For each q = p, . . . ,1 + m + n, the support of V contains exactly one triplet (i,j,k) with 
% + j + k = q. 

2. Those triplets are well-ordered. (The partial order given in (pQ) is a total order on these triplets.) 
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3. All other triplets in the support have entries which are index-wise smaller than or equal to those 
of the unique triplet (io,jo, ko) with iq + jo + ko = p. 

For example, the only vertex "well-ordered starting at level 3" is the well-ordered vertex V. Slightly 
less trivially, it is also the unique vertex "well-ordered starting at level 4" . On the other extreme, all 
vertices that contain (l,m,n) as a support triplet are well-ordered starting at level l + m + n. Observe 
that from any vertex of T x>yjZ we can move, by a single pivot edge in the sense of the simplex method, 
to another vertex containing any prescribed entry (i,j,k) to be non-zero. In particular, we can move 
to a vertex that has (I, m, n) in its support. So, we can assume from the beginning that (I, m, n) is in 
the support of our initial vertex V , and will add one to the count of edges traversed to arrive to V . 

Lemma 4.4 IfV is a vertex ofT X)V ^ z that is well-ordered starting at level p € {5, . . . , l + m + n}, then 
there is a path of at most 2(p — 4) edges of T X) y^ z that leads from V to a vertex that is well-ordered 
starting at level p — 1 . 

Proof. Let (io, jo, ko) be the unique triplet in the support of V with io + jo + &o = P- We first observe 
that there is no loss of generality in assuming that p = l + m + n (that is, {jo, jo, ko) = (I, m, n)). This 
is because the vertices of T x>y>z that are well-ordered starting at level p and agree with V in all the 
triplets with sum of indices greater than or equal to p are the vertices of a non-degenerate io x jo x ko 
axial transportation polytope, obtained as in the proof of Lemma l4.1i 

So, from now on we assume that V is well-ordered starting at level p = l + m + n. Let Si be the 
set of support triplets in V, other than (l,m,n), that have first index equal to I. Similarly, let S2 and 
S3 be the sets of support triplets that have, respectively, second and third indices equal to m and n. 

Our goal is to modify V until one of Si, S2 or S3 becomes empty, but always keeping the triplet 
(l,m,n) in the support. Once this is done, a single pivot step can be used to obtain a vertex that 
is well-ordered starting at level p — 1 as follows: Without loss of generality assume that Si is empty 
(the cases when S2 or S3 are empty are treated identically). In particular, neither (l,m — l,n) nor 
(/, m,n — l) are in the support. If (I — 1, m, n) is in the support then our vertex is already well-ordered 
starting at level p — 1. If not, we do the pivot step that inserts (I — l,m,n). This pivot step cannot 
remove (l,m,n) or insert (l,m — l,n) or (l,m,n — 1) in the support. (The (l,m,n) coordinate is 
not removed from the support since the entry remains constant in this pivot. The (l,m — l,n) and 
(l,m,n — 1) coordinates remain zero because only non-zero entries of V and the entry (I — l,m,n) 
change in the pivot). This pivot produces a vertex well-ordered starting at level p — 1. Figure H] gives 
a picture for this case. 

Given a vertex V well-ordered at level p, we specify a sequence of pivots in the graph of T X)y>z to 
a vertex V' such that one of Si, S2 or S3 is empty for V. Lemma 14.51 below shows how to get such a 
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Figure 4: A well-ordered vertex starting at level p — 1. 

V in a number of steps bounded above by 

2|Si U S 2 U S 3 | - 3 < 2(p - 3) - 3 = 2p - 9. 

In one more step, that is, at most 2p — 8, we get to a vertex that is well-ordered starting at level p — 1. 
This completes the proof of our lemma. H 

For Lemma 14.51 let us introduce the following notation: 

R l :=S 1 \(S 2 US 3 ), R2 :=S 2 \(S l US 3 ), R 3 := S 3 \ (Si U S 2 ), 

#12 : = £l H 52, -Rl3 := 5i Pi 53, i?23 : = 5 2 fl 53- 

That is, Ri consists of the elements of S1US2US3 that belong only to Si, and Rij of those that belong 
to Si and Sj. Observe that, by definition, no element belongs to the three Si's, so that S± U £2 U S3 
is the disjoint union of these six subsets. 

Lemma 4.5 With the above notation and the conditions of the proof of the previous lemma, suppose 
that no Si is empty. Then: 

1. If both Ri and Rjk are non-empty, with {i,j,k} = {1,2,3}, then there is a single pivot step that 
decreases \Si U S2 U S3I. 

2. If the three Rij 's are non-empty, then there is a sequence of two pivot steps that decreases 

\s 1 us 2 us 3 \. 

3. If the three Ri 's are non-empty, then there is a sequence of two pivot steps that decreases \S\ U 

s 2 us 3 \. 

4- If none of the above happens, then S\ U S 2 U S3 is contained in one of the Si 's, say S\. Then, 
there is a sequence of \S± \ — 1 pivot steps that makes S± U S2 U S3 empty. 
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All in all, there is a sequence of at most 2\S\ U S2 U Ss\ — 3 pivot steps that makes some Si empty. 

Proof. Let us first show how the conclusion is obtained. We argue by induction on \Si U S2 U S$\, the 
base case being \S\ U S2 U S^\ = 2, which is the minimum to have no Si empty. The base case implies 
we are in the situation of either part (1) or part (4), and a single pivot step makes an Si empty. 

If \S\ U 52 U S3I > 2 and one of the conditions (1), (2) or (3) holds, then we do the step or the 
two pivot steps mentioned there and apply induction. If none of these three conditions hold then 
it is easy to see that (4) must hold. (Remember that we are assuming that no Si is empty, and 
Si = Ri U Rij U Rik)- Part (4) guarantees we have a sequence of |5*i| — 1 < 2|5i U S2 U ^3 1 — 3 pivot 
steps that makes an Si empty. 

So, let us prove each of the four items in the lemma. Let a denote the entry (l,m,n). 

1. Suppose without loss of generality that R\2 and R3 are not empty. Let (3 G R12 and 7 € R3. The 
reader may find it useful to follow our proof using Figure [5] which depicts the situation. Observe 




Figure 5: A layout of entries a, (3, 7 and 5. 



that a = (I, m, n) is the index-wise maximum of [3 and 7. Let 5 be the index-wise minimum 
of them. First observe that 5 is not in the support of vertex V. Otherwise, we could add 
±| minjao,, ap, o 7 , as, } (e a + e$ — ep — e 7 ) to V and stay in T XtVtZ . Hence, V would be a convex 
combination of two other points from T X mz (and thus not a vertex), parallel to the direction of 
e-a + e<5 — ep — e 7 (here eij,u denotes the basis unit vector in the direction of the variable aj^fe). 

Next, consider V' = V + minjag, a 7 } (e a + e$ — ep — e 7 ). Observe that V has different support 
than V since either (3 or 7 has been removed (not both, because ap ^ a 7 by non-degeneracy). 
Also, since V cannot have support strictly contained in that of V, 5 must have been added. 
That is, the supports of the vertices V and V differ in the deletion and insertion of a single 
element, which means they are adjacent in the graph of the polytope T x>ytZ . As desired, when 
going from V to V' the cardinality of S\ U S2 U S3 is decreased by one. 
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Example 4.6 (Example 14.31 continued) Consider the vertex V with non-zero coordinates 

«(1,1,2) = 28 ; a (2,l,2) = I 2 ; a (2,2,3) = 6 , 1(1,3,2) = 2 , °(1,3,1) = 82 > a (3,3,2) = 2 > and ^3^3) = 28. 

In £/ws example, a = (3,3,3), /3 = (3,3,2), 7 = (2,2,3), and 5 = (2,2,2). After clearing ag, 
we arrive at the vertex V with non-zero coordinates a'^ 1 2) = 28 > a> (2 1 2) = ^ 2 > a (2 23) = ^> 

a (l,3,2) = 2 ' a (l,3,l) = 82 ' a (2,2,2) = 2 ' anrf a (3,3,3) = 28> 

2. Suppose now that none of the Rij's is empty, and let G i? 13 , 7 € -R23 and <5' G i?i2- we apply 
the same pivot as in case one, which makes 5, the coordinate-wise minimum of f3 and 7, enter 
the support. This pivot does not decrease \S± U 5*2 U S3I, but it leads to a situation where we 
have S G R3 and 5' G i?i2- Hence, we can apply part one and decrease |Si U S2 U S3I with a 
second step. 

4. Let us prove now part (4) and leave (3), which is more complicated, for the end. Observe that if 
Si U S2 U S3 = Si but S2 and S 3 are not empty, then necessarily R12 and i?i3 are both non-empty. 
While this holds, we can do the same pivot steps as before with a (3 G -R12 and a 7 € i?i 3 . Each 
step decreases by one the cardinality of -R12 U -R13, increasing the cardinality of R\. The process 
finishes when R\2 (hence S2) or i?i 3 (hence S3) becomes empty, which happens, in the worst 
case, in |Si| — 1 steps. 

3. Finally, consider the case where the three Ri's are non-empty. Let = (l,ji,k\) G R\, 7 = 
(22,^,^2) G R2 and 5 = (i3,j3,re) G R3, and, as before, a = (l,m,n). Figure [6] depicts the 
situation to help the reader with following our proof. Let ei and €2 be two triplets of indices 
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Figure 6: Possible layouts of the entries a, (5, 7, 5, e\ and €2- 

with the property that {a, ei, £2} and 7, 5} use exactly the same three first indices, the same 
three second indices, and the same three third indices. For example, let us make the choice 
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ei = (*2)ii,fci) and €2 = (13,^3,^2) as in the left part of Figure El By non-degeneracy, the 
smallest value among ap, a 7 and as at V is unique. We assume without loss of generality that 
the smallest among them is ap. Let W be the point of T Xj y jZ obtained by changing the following 
six coordinates: 

a' a = a a + ap, alp = ap — ap = 0, a'^ = a 7 - ap, a' s = as - ap, a' ei = a €1 + ap, a' e2 = a 62 + ap. 

It may occur that ei = 62, as shown in the right side of Figure EJ Then we do the same pivot 
except we increase the corresponding entry a ei = a e2 twice as much. 

Observe that one of e% or €2 may already be in the support of V, but not both: Otherwise W 
would have support strictly contained in that of V, which is impossible because V is a vertex 
and has minimal support. If one of ei or 62 were already in the support of V, or if ei = 62, then 
W is a vertex and we take V = W. As in the first case, V is obtained from V by traversing a 
single edge and has one less support element in Si U 52 U S3 than V: None of e\ or €2 can have 
a common entry with a, since none of f3, 7 and 5 has two common entries with a. 

However, if e\ and €2 are different and none of them was in V, then W has one-too-many 
elements in its support to be a vertex, which means it is in the relative interior of an edge 
E and L = VW is not an edge. See Figure [71 Moreover, both E and VW lie in a two- 
dimensional face F. This is so because every support containing the support of a vertex defines 
a face of dimension equal the excess of elements it has. In our case, F is the face with support 
support(y) U support(W) = support(y) U {ei, 62}- 



E 




E 



Figure 7: The octagon containing segment VW arising when entries e\ and €2 are different and none 
of them was in the support of V 

We now look more closely at the structure of F. Each edge H of F is the intersection of F with 
a facet of our transportation polytope. That is, there is a unique variable 77 that is constantly 
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zero along H but not zero as we move on F in other directions. For example, since e\ and €2 
are zero at V but not constant on F (they increase along L), V is the common end of the edges 
defined by ei and €2- 

Our goal is to show that there is a vertex V' of F at distance at most two from V and incident to 
the edge defined by one of the variables j3, 7 and 5. At such a vertex V' we will have decreased 
[Si U S2 U S3 1 by one, as claimed. The key remark is that there are at most two edges of F not 
produced by one of the variables a, (3, 7, 5, e%, and €2- every variable r\ other than those six is 
constant along L, so it either produces an edge parallel to L or no edge at all. In particular, F 
is at most an octagon, as in Figure Now: 

• If F has five or less edges, then every vertex of F is at distance one or two from V. Take as 
V' either end of the end-point W of L. This works because at W one of /3, 7 or 5 is zero, 
by construction. 

• If F has six or more edges, then the two vertices V and V" of F at distance two from V 
are at distance at least two from each other. So, together they are incident to four different 
edges, none of which is the edge of e\ or 62- (Remember that the edges of ei and €2 are 
incident to V). At least one of these four edges is defined by 0, 7, or 5, because there are 
(at most) three other possible edges: the one of a and two parallel to L. 

H 



Example 4.7 To make ideas completely clear, using the same polytope T x ,y,z as in Example \4-3\ we 
consider its vertex V with non-zero coordinates ^(1,1,3) = 25, ^(3,1,1) = 15, 1^3, 2 ,i) = 6, ^(1,3,1) = 61, 
v (i,3,2) = 26, ^(2,3,2) = 18 and ^(3,3,3) = 9 - Here , « = ( 3 > 3 > 3 ); I 3 = (3, 2, 1), 7 = (2, 3, 2), 5 = (1, 1,3), 
ei = (2,2,1), and 62 = (1,1,2). The triplet (3 is not in the support of W, and W has non-zero 
coordinates ^(1,1,3) = 19, ^(1,1,2) = 6, u> (3>lil ) = 15, ^(2,2,1) = 6, W(i Al) = 61, ^(1,3,2) = 26, ^(2,3,2) = 
12, and ^(3,3,3) = 15- 

The vertices ofT X)VyZ with support contained in support(V) U {e\,€2} form the A-gon F = VBCD 
where B is the vertex with non-zero coordinates 

k(i,i,3) = 12,fe(i i i i 2) = 26,6( 3ill ) = 2,6(32,1) = 6,6( 3i 3 3) = 22,6( 2i3i 2) = 18, 6(1,3,1) = 74, 
C is the vertex with non-zero coordinates 

c (l,l,3) = 22, £(3,1,1) = 18, C( 2 ,2,l) = 6, C(3,3,3) = 12, C(i,3,2) = 32, C( 2 ,3,2) = 12, C(l,3,l) = 58 
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and D is the vertex with non-zero coordinates 



d(l,l,3) = 6,(^(1,1,2) 



32,d( 3jljl ) = 2,d( 2i2 ,i) = 6,d( 3j3i 3) = 28, ^(2,3,2) = 12,d(i,3,i) = 74. 



Note W is in the edge E = CD. We let V' = D, the endpoint of E closer to V. Thus, we use one 
edge to go from V to V . 

Proof of Theorem \1.1\ Starting with any vertex "well-ordered starting at level p = I + m + n" (which 
can be reached in a single step) we use Lemma 14.41 to decrease one unit by one the level at which our 
vertex starts to be well-ordered until we reach the unique well-ordered vertex V . Thus, the number 
of steps needed to go from an arbitrary vertex V to V is at most 



To go from one arbitrary vertex to another, twice as many steps suffice. H 

Remark: The whole proof can be generalized to arbitrary axial d-way tables, instead of d = 3, 
without much effort. Everything in Lemma [4.11 goes through without change, as well as the definition 
of "well-ordered starting at level p" . In the other arguments, the first change is that we have d sets 
Si, . . . , Sd instead of just three. In particular, in the proof of Lemma 14.51 the worst case will be that 
of d — 1 different e's, which gives a face of dimension d — 1. Hence, the bound given in the statement 
of Lemma 14.41 can be substituted to the maximum diameter of a simple polytope of dimension d — 1 
with at most p facets. This still yields a polynomial bound for any fixed value of d. We leave the 
details for the interested reader. 
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